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Abstract 

We investigate the scattering of phase oscillation of Bose-Einstein 
Condensate by a 'draining of bathtub' type fluid motion. We derive 
a relation between the reflection and transmission coefficients which 
exhibits existence of analog of 'Superradiance effect' in BEC vortex 
with sink. 

1 Introduction 

It has been a hundred years since Einstein's discovery of general relativity 
which has been experimentally verified to a high degree of accuracy. However, 
on the semiclassical front the absence of direct experimental observation of 
the phenomenon of particle emission predicted by Hawking [1 in 1975, still 
stands out as a problem. In 1981 though, Unruh Ej showed that if a 
fluid is barotropic and inviscid, and the flow of the fluid is irrotational, the 
equation of motion that fluctuation of the velocity potential of acoustic dis- 
turbance obeys, is identical to that of a minimally coupled massless scalar 
field propagating in an effective curved space-time Lorentzian geometry. One 
can simulate artificial black hole 0] in this system by trapping the sound 
wave in a compact region using radial velocity of the fluid flow. Introduc- 
tion of a rotational component [H 03 IE] to that radial velocity will make the 
artificial black hole an analog of rotating black hole in physical gravitational 
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system. He also pointed out that there is possibility of experimental obser- 
vation of the acoustic analog of Hawking radiation from regions of flow of 
inviscid and barotropic fluids behaving as outer trapped surfaces ('acoustic 
event horizons'). This model is sufficiently rich to enable probing of almost 
all the kinematic aspects of general relativity but not its dynamics, because 
dynamics of a fluid system is governed by Euler equation and continuity 
equation but dynamics of a gravitational system is governed by Einstein's 
equations. 

In our earlier work [7J |H], using the analogy between a shrinking fluid 
vortex ('drain bathtub') modelled as a (2+1) dimensional fluid flow with a 
sink at the origin and a rotating (2+1) dimensional black hole with an ergo- 
sphere, it has been shown that a scalar sound wave is reflected from such a 
vortex with an amplification for a specific range of frequencies of the incident 
wave, depending on the angular velocity of rotation of the vortex. This is 
analogous to superradiant scattering [31101111] which occurs in the case of ro- 
tating black holes in asymptotically flat spacetime. A more detailed analysis 
this phenomenon in the low frequency range has also been done in the next 
paper The reflection coefficient exhibit its frequency dependence which 
is consistent with the results in jH] . We have also discussed the possibility of 
observation of this phenomenon especially for inviscid fluids like He-II, where 
vortices with quantized angular momentum may occur. 

Although we have demonstrated the analog of superradiance in the su- 
perfluid helium, BEC is a far cleaner system which can be manipulated to 
produce the kind of metric we have (DB) in a much easier manner. That is 
why BEC seems to be better suited for this purpose. Recent experimental 
investigation of Bose-Einstein condensation ^21 EH E3] has further improved 
the prospects of experimental realization of the analog of kinematical effects 
of gravity. 

In this paper we investigate the possibility of the acoustic analog of super- 
radiance scattering from hydrodynamic vortex with a sink in Bose-Einstein 
condensate. Linear acoustic perturbations in such a system are shown to scat- 
ter from the ergoregion with an enhancement in amplitude, for a restricted 
range of frequencies of the incoming wave. 

The paper is organised follows. First, in section(2), we briefly discuss 
how effective Lorentzian geometry emerges from the propagation of phase 
oscillation in BEC. In section(3) we briefly describe the properties of effec- 
tive geometry that emerges from propagation of phase oscillations through 
BEC. In section(4) we demonstrate analog of superradiance effect in BEC 
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for position dependent background configuration. We conclude with some 
remarks in section(5). 



2 Effective Geometry in BEC 

A Bose-Einstein condensate |15| ITB] is a collection of indistinguishable bosons 
such that all of them are in same single-particle quantum state 0(f). If there 
is no interaction between the bosons, then the second quantized Hamiltonian 
for N such bosons is given by, 



N 



and the corresponding wave function of this system is, 



$(fi,r 2 r N ) = cj)(ri) <f>(f 2 ), ,^>(r N ) (2) 

Even there are very week interaction between the atoms, at 0°k almost all the 
atoms in the condensate remains in the same single-particle state. If the number of 
boson is large enough then the time evolution of such a condensate is determined 
by Gross-Pitaevskii equation, 

m*=(-^ + v^ + ^\*A*, (3) 

V 2m m I 

where m is the mass of the atoms, a is the scattering length, and we normalise 
to the total number of atoms J d 3 x|<&(x, t)\ 2 = N. 

Expressing $(f, t) in terms of its amplitude \J p(f, t) and phase VP(r, t) , 
the above wave equation© is converted into two first order partial differential 
equation, 



dt* = -— (Vtf.Vtf)-V p+- -V 2 ^ 

2m m 2m yjp 

dtp = - — (Vp.V^ + pV 2 ^) (4) 
m 

Now linearising p and ^ around there background values po and \&o respec- 
tively, 
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9t*i = — v^o-v^i p + — tt^v 2 -^L -pi— 

m m 2m {2^ \\fPo J 2 pz J 

d*pi = V(poVfi) V(pi Wo) 

m m 

Here p\ and are the perturbed values of p and respectively. 
If we neglect the quantum pressure term, then the combination of above equa- 
tions© results the following covariant differential equation, 



where, 



1 I _1 ! " V o 



The speed of sound(c) and the background velocity(uo) are defined as, 

c(r) = —JA7rap (r) vi = — d;# (7) 

m v m 

Therefore the propagation of the phase oscillations in Bose-Einstein condensate 
is same as the propagation of a minimally coupled massless scalar field in an 
effective Lorentzian geometry which is determined by the background velocity 
and background density of the condensate and also by the speed of sound in the 
condensate [T71 UK) IP3] . 

Since in this paper, we investigate the possibility of the acoustic analog of 
superradiance (a phenomenon that we call 'superresonance'), i.e., the amplification 
of a sound wave by reflection from the ergo-region of a rotating acoustic black 
hole, we choose the so-called 'draining bathtub' type of fluid flow [3], which is 
basically a (2+1) dimensional flow with a sink at the origin. A two surface in this 
flow, on which the fluid velocity is everywhere pointing towards the sink, and the 
radial velocity component exceeds the local sound velocity everywhere, behaves 
as an outer trapped surface in this 'acoustic' spacetime, and is identified with the 
(future) event horizon of the black hole analog. Thus, the velocity potential for 
the flow has the form 5 (in polar coordinates on the plane), 

ip(r,4>) = A logr + B <j> , (8) 
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where, A and B are real constants. This leads to the velocity profile, 



A. B 7 

v = r H <p . (9) 

r r 

For the velocity potential given in equation©, the analog black hole metric is 
(2+1) dimensional with Lorentzian signature, and is given by, 



ds 2 = ^ 



A 2 + B r 



dt* 



2A , , 
— ar at 

r 



2Bi 



> dt + dr 2 + r 2 i 



(10) 

where c is the velocity of sound. In our earlier work we had neglected the 
conformal factor in the effective metric because the back ground density and hence 
the background pressure of fluid were assumed to be constant. In this paper 
we have relaxed this condition by considering position dependence of the density 
profile, which leads to the position dependence of the speed of sound also. 

Recently Slatyer and Savage [20] have suggested the following density profile 
|16j which could be realized(at least nearly) in future BEC experiments, 



Po(r) = Po 



ro) 



(11) 



(r - r ) 2 + 2cr 2 

where a is the spatial scale of variance, the distance over which the wave 
function of the Bose-Einstein Condensate tends to it bulk values, when the wave 
function is subjected to a localized perturbation. For analytical simplicity we have 
considered ro = case in this paper. The sound speed© corresponding to this 
density profile is, 



c(r) 



- 2 + 2a 2 



where, 



— y/4wap 
m 



(12) 



(13) 



Considering the above velocity profile© and new density profile yields new 
locations of analog of event horizon and ergosphere. As for the Kerr black hole 
|21j in general relativity, the radius of the ergosphere is given by the vanishing of 
goo, i-e., at 



(A 2 + B 2 
2 r 2 



1 + 1 + 



8cla 2 
A 2 + B 2 



(14) 
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The metric has a (coordinate) singularity at, 



rh 



A 2 

2 



1+ 1 + 



8c^ 2 
A 2 



(15) 



, which signifies the horizon, i.e., the boundary of the outer trapped surface. 

It should be noticed that both the radius of the horizon and the ergosphere are 
larger than what we found in our earlier work assuming constant density profile. 
This result is expected because in this case the speed of sound is decreasing towards 
the centre of the vortex and hence it becomes lesser than fluid velocity at larger 
value of radial coordinate. 



3 Superresonance Scattering 

From the components of the draining vortex metric it is clear that the (2+1)- 
dimensional curved spacetime possesses isometries that correspond to time trans- 
lations and rotations on the plane. The solution of the massless Klein Gordon 
equation© can therefore be written as, 

#(i, r, (f>) = Re{R(r) e' 1 w * j m (16) 

where u and m are real and positive. In order to make \P(i, r, (j>) single valued 

, m should take integer values. 

Then the radial function R(r) satisfies, 

f^M + P l(r) ^) +gi(r)fl(r) = 0, (17) 

where, 



P x {r) 



r 2 + 2a 2 



Coor {r 4 (r 2 h + 2a 2 



- 2 + 2<j 2 )} 



P2{r) 



R: 



'2U- 
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Coo r h 



2iWr 2 rlxlr 2 h + 2a 2 + 



(r 2 + 2a 2 ) (r 2 + 6a 2 ) 
(r 2 + 2a 2 ) 2 



Qi(r) 



r 2 + 2a 2 



c^r 2 



(r 2 — r%) (r 2 r\ + 2 (r 2 + r 2 ) a 2 ) 



Q 2 (r) 
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g 2 (r) = 2iB Coo mr 2 h (r 2 + 2<r 2 ) + r 4 ^ + 2<r 2 [ Coo 2 m 2 -W 2 (r 2 + 2<r 2 )} 
and, 

(18) 

Now in order to eliminate the imaginary part from the P\{f) and Q\{r) we 
made the following substitution, 

R(r) = e ^/iM e -imMr) Q{r) (19) 



r 2 +2a 2 



fl( r ) r. f 4 i c 2 9 i o 4\ fs(. r ) 

2 Coo (r£ + 6 r£ a z + 8 a 4 ) 



A(r) = 2^nE(fj + 4 ff ») A(r) 



AM = (rJ + 2,7 2 ) 2 log(r 2 ->f) -4<j 4 log{r 2 rJ + 2 (r 2 + rj) a 2 } 



am - + 2 . 2 ) io g (^i) + 2 . 2 log p^ +2 '; +r s- 2 ) 

In terms of new radial function G(r) the radial eauation ()17() takes the following 
form, 

d -^r- + Sl{r) + S2{r) G{r) = (20) 

where, 



3 r 4 + r\ 2 r 2 r r\ 

r (r 4 — r£) r 2 + 2 a 2 (r 2 + r^) (r 2 r 2 + 2 (r 2 + r 2 .) a 2 ) 



a r \ r 2 (r 2 +2a 2 ) 

Coo ( r r h ) V rf l + 2r z o l + 2rf l a z ) 
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S 3 ( r ) = W 2 r* (r 2 + 2a 2 ) (r 2 h + 2 a 2 ) + & m 2 [r\ (r 2 + 2 a 2 ) - r 4 (r 2 + 2a 2 )} 
Now we introduce tortoise coordinate r* through the equation, 

d ( A 2 \ d 



dr* \ r 2 c 2 I dr 

which implies that, 



(21) 



2r h a 3 ,( V2r h a \ r h (r 2 + 2 a 2 ) tanh 



cot 



± + a 2 (r 2 +4a 2 ) \r^r 2 + 2a 2 ) r 2 + 4a 2 

This tortoise coordinate spans the entire real line as opposed to r which spans 
only the half-line. The horizon at r = maps to r* — > — oo , while r — > 
oo corresponds to r* — > +oo . 

Let us now define a new radial function G(r) as, 



G(r) = Vr2 \ 2a2 H(r) (22) 

T 2 

Substituting this in equation© we observe that H{r) satisfies the differential 
equation, 

, . d 2 H(r) , . dH(r) , x , , 
P (r) , I' + Q(r) + V(r) H{r) = 23 

where, 



P(r) 



j (r 2 + 2a 2 ) ] 2 
-« (r 2 + 2a 2 ) f 



_ 2r 4 (r 2 + 4a 2 ) f _ r{ (r 2 + 2a 2 ) 
^ [ ' (r 2 +2a 2 ) \ r 4 (r 2 +2a 2 ) 



v(r) = (r 2 + 2* 2 ) w2 (r 2 -r 2 ) (r 2 r 2 + 2 r 2 a 2 + 2 r 2 a 2 ) 
c^r 2 (r 2 + 2a 2 ) 2 {r 2 + 2a 2 ) 2 

and, 
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51 r\ a 4 42r 4 a 6 a 2 {35r 4 + 2a 2 (3 + 4m 2 ) (r 2 + 2a 2 )} 
^v) — ~s 



+ 



r 8 j,10 2 r§ 

1 - 4m 2 ) (r 2 + 2a 2 ) 5r 4 - a 2 (5 - 4m 2 ) (r 2 + 2a 2 ) 



4 t-2 4r^ 



Now in terms of tortoise coordinate one obtains the modified differential equa- 
tion as, 

d 2 H(r*) 

—^1 + V(r)H(r*) = (24) 

We analyse this differential equation in two distinct radial regions, viz., near 
the sonic horizon, i.e., at r* — ► — oo and at asymptopia, i.e., at r* — > +oo . 

Our choice of coordinates and substitutions are with a view to transform differ- 
ential equation (j!7|) into a forml|24|) which is simple enough to identify ingoing and 
outgoing modes both in asymptotic and near horizon regions. In addition we have 
managed to eliminate the first derivative term from the differential equation which 
leads to constancy of Wronskian of the differential equation. This will facilitate 
the calculations of the transmission and reflection coefficients. Since the second 
term of V(r) in eauation ()24j) tends to zero both in asymptotic region and near the 
horizon, the first term in V(r) is the only important term. 

In the asymptotic region, the above differential eauation H24j) can be written 
approximately as, 

d 2 H(r*) uj 2 



oc 



This can be solved trivially, 

H(r*) = n wm exp (i ^- r*^j + exp (-i ^- r*^j (26) 

The first term in equation (J^HJ) corresponds to reflected wave and the second 
term to the incident wave, so that R is the reflection coefficient in the sense of po- 
tential scattering. It is not difficult to calculate the Wronskian of the solutions Q26JI : 
one obtains, 

W(+oo) = 2 i — (l - |7^ m | 2 ) . (27) 



c 



>3C 



Similarly , near the horizon the above differential equation can be written 
approximately as, 
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d 2 H(r*) 

1 + 



dr 



*2 



(r 2 h + 2a 2 ) (uj-mQ H f 



H{r*) = (28) 



where £Ih is the angular velocity of the sonic horizon. We impose the physical 
boundary condition that of the two solutions of this equation, only the ingoing one 
is physical, so that one has, 



H(r*) = % m exp { -i ,1 M±^i r * } (29 ) 



Once again, it is easy to calculate the Wronskian of this solution; one obtains, 

) 



W(-oo) = 2 i ( ^ \% m \ 2 . (30) 



Ch 



Now using these approximate solutions (|2H|29|) of the above differential equa- 
tion together with their complex conjugates and recalling the fact two linearly 
independent solutions of this differential equation l)24|l must lead to a constant 
Wronskian, it is easy to show that, 



where, 



1 - |7^ m | 2 = ( 1 - ^) °— \% m \ 2 (31) 
1 uj J c h 



c " = (it+i?) ' c ~ (32) 



Here 1Zumi and 7^ m are the amplitudes of the reflection and transmission 
coefficients of the scattered wave, respectively. It is obvious from equation (|31|) 
that, for frequencies in the range < uj < m £Ih, the reflection coefficient has 
a magnitude larger than unity. This is similar to the amplification relation that 
emerges in our earlier analysis of superresonance |S] . The only difference here 
is the appearance of a factor Cqo/c^ . This is due to the fact that speed of sound 
in Bose-Einstein condensate depends on radial coordinate. 



4 Conclusion 

In this paper we have shown that when acoustic disturbances propagate through 
a Bose-Einstein condensate, they get scattered by the flow of the condensate and 
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in the particular frequency range acoustic disturbances are reflected with an am- 
plitude greater than the amplitude of the incident wave. The frequency range is 
determined by angular momentum of the incident wave and angular velocity of the 
horizon. The extra energy that reflected waves carry is taken from the rotational 
energy of the condensate. As a result vortex motion is slowed down and eventu- 
ally stops when all the rotational energy of the condensate is extracted out of it. 
At other frequencies acoustic disturbances transmit through the condensate and 
ultimately get trapped inside the Dumb hole. This is analogous to the superradi- 
ance effect that occurs exclusively in the case of rotating black holes in physical 
gravitational systems. Because of the discreteness of the parameter B of rotational 
component of the velocity profile (JSJ)(as evident in BEC), the angular momentum 
of the acoustic black hole will also be proportional to an integer and hence causes 
the energy flux to change discretely making observation of the phenomenon easier 
than that in actual black holes. This also manifests itself in the reflection coeffi- 
cient whose spectrum has equally spaced peaks of varying strengths all of which 
are multiples of the minimum strength. These interesting aspects are worth being 
studied quantitatively in greater detail. 

Unlike our earlier work we have considered, following the suggestion of Slatyer 
and Savage [2111 > a particular type of spatial dependence of the background density 
of the condensate. Presence of spatial scale of variation of the condensate wave 
function has nontrivial effect on the artificial Lorentzian geometry that emerges 
from BEC. Due to spatial dependence of the background density of the fluid, speed 
of sound in such a system varies when it propagates though the fluid. This spatial 
dependence modifies the relation between reflection and transmission coefficients 
of the scattered waves by a factor Cqo/c^, but does not affect the cutoff frequency 
of the superradiance effect. In the limit a tends to zero the background density 
and hence the speed of sound in the fluid become constant. As a check we note 
that in this limit we recover the results of our earlier works . 
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